Anisotropy and non-universality in scaling laws of the large scale energy spectrum in 

rotating turbulence 
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Rapidly rotating turbulent flow is characterized by the emergence of columnar structures that 
are representative of quasi-two dimensional behavior of the flow. R is known that when energy is 
injected into the fluid at an intermediate scale Lf, it cascades towards smaller as well as larger 
scales. In this paper we analyze the flow in the inverse cascade range at a small but fixed Rossby 
number, IZof ~ 0.05. Several numerical simulations with helical and non-helical forcing functions 
are considered in periodic boxes with unit aspect ratio. In order to resolve the inverse cascade 
range with reasonably large Reynolds number, the analysis is based on large eddy simulations which 
include the effect of helicity on eddy viscosity and eddy noise. Thus, we model the small scales 
and resolve explicitly the large scales. We show that the large-scale energy spectrum has at least 
two solutions: one that is consistent with Kolmogorov-Kraichnan-Batchelor-Leith phenomenology 
for the inverse cascade of energy in two-dimensional (2D) turbulence with a ~ fc^ scaling, and 
the other that corresponds to a steeper ~ kj 3 spectrum in which the three-dimensional (3D) modes 
release a substantial fraction of their energy per unit time to the 2D modes. The spectrum that 
emerges depends on the anisotropy of the forcing function, the former solution prevailing for forcings 
in which more energy is injected into the 2D modes while the latter prevails for isotropic forcing. 
In the case of anisotropic forcing, whence the energy goes from the 2D to the 3D modes at low 
wavenumbers, large-scale shear is created resulting in another time scale t 3 h, associated with shear, 
thereby producing a ~ fc _1 spectrum for the total energy with the 2D modes still following a ~ k± 
scaling. 

PACS numbers: 47.32.Ef, 47.27.Gs, 47.27.Jv 



I. INTRODUCTION 

The emergence of anisotropy in rotating flows has been 
studied extensively since the seminal work of Taylor and 
Proudman [JH5]. More recently, it has been observed 
both experimentally [U |5] and numerically [6] for the 
velocity field, as well as for passive scalar fluctuations 
using a reduced model [7J [5] and direct numerical sim- 
ulations (DNS) ^9 . While linear theory like the Taylor- 
Proudman theorem (see, e.g., |10j ) explains the existence 
of the columnar structures in the laminar regime, more 
recent theories account for the nonlinearities in the flow 
and shed light on the mechanism of two dimensionaliza- 
tion of the flow from an initial isotropic state [TTUI3] . 



An interesting feature of such flows is the trans- 
fer of energy to large scales on application of exter- 
nal forcing [14TU6| . However, there has been little 
consensus about the scaling of the energy spectrum 
at large scales. The numerical simulation of [TS] re- 
spectrum whereas [141 116] show a 



ports a steep ~ k ± ' 
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more conventional ~ k ± scaling that is reminiscent of 
two-dimensional (2D) Kolmogorov-Kraichnan-Batchelor- 
Leith (KKBL) phenomenology for an inverse cascade of 
energy [33] (also see [H]). In [IS], a model is used to 



show that a ~ k ± spectrum for the 2D modes (also 
called "slow modes" ) results when triadic interactions be- 
tween the 2D and the 3D ("fast") modes are discounted 
for artificially, but a ~ fej spectrum is observed when 
all interactions are accounted for. Besides, a ~ kj 3 law 
for the horizontal kinetic energy spectra is also observed 
in a rapidly rotating Rayleigh-Benard convection using 
a reduced model [TS]. It must be pointed out that the 
conserved quantities (for an inviscid fluid) in the case of 
two dimensional flows, viz., energy and enstrophy, are 
different from the three dimensional case, where energy 
and helicity (the correlation between velocity and vortic- 
ity) are conserved. Therefore, the physical mechanism 
leading to an inverse cascade of energy in the three di- 
mensional case does not follow immediately from its two 
dimensional counterpart. Besides, the nature of forc- 
ing, i.e. spectrally isotropic vs. anisotropic (equivalently 
3D vs. 2D), with three or two spatial components (3C 
vs. 2C), and the aspect ratio of the computational box 
may play significant roles in the dynamics of the flow. 

In this paper, we revisit the issue of inverse cascade of 
energy in rotating flows within a specific framework, viz., 
fixed Rossby number and unit aspect ratio of the compu- 
tational box. However, we explore different forcing func- 
tions to consider the effects of spectral anisotropy and 



of helicity in the inverse cascade range. We present re- 
sults from numerical simulations that use a subgrid scale 
model developed in [5D] ; this model was validated against 
DNS of rotating flows in [2"TH2"5] . We show that the two 

observed spectra, viz., ~ kj_ and ~ fej can arise in 
full simulations (simulations that resolve all triadic inter- 
actions and account for coupling between the 2D and the 
3D modes). When the forcing is isotropic, energy goes 
from the 3D to the 2D modes and a ~ kj 3 spectrum 
results for the energy in the slow modes. When more en- 
ergy is pumped into the 2D modes, less energy goes from 
the 3D to the 2D modes and a ~ fc~ ' spectrum is seen 
for the slow modes. 

The kinematics of the nonlinear advection term also 
changes significantly as the spectra of slow modes change 



from 



fcj/ to 
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We study the velocity gra- 



dient tensor in all simulations and compute the largest 
eigenvalue of the rate of strain tensor. For the case of 
anisotropic forcing when the flux of energy between the 
3D and the 2D modes reverses and energy at large scales 
goes from the 2D to the 3D modes, a significant amount 
of shear is created at large scales. This introduces a new 
shear timescale r s h that is independent of wavenumber. 
As a result, the spectrum for the total energy approaches 
a ~ fc _1 power law. 

The remainder of this paper is organized as follows. In 
Sec. In] we discuss previous results, introduce equations 
and notations used in the rest of the paper, and derive 
equations to study the coupling between modes and the 
energy transfer between scales. In Sec. |III| we present 
the LES model used in the numerical simulations and 
describe all the runs as well as the different spectra used 
to characterize scaling laws in the inverse cascade range. 
Finally, in Sec. |IV| wc present and discuss the numerical 
results, while in Sec. Wj we conclude with brief remarks 
and pointers to some open questions. 



II. INERTIAL WAVES AND ENERGY 
TRANSFER TO THE SLOW MANIFOLD 

A. Equations 

The non-dimcnsionalized incompressible Navier-Stokes 
equations with global rotation, $1 = £lz, are as follows: 



<9 t u+(u-V)u+ — ixu=-VP+ — V 2 u + f, (1) 
TZo TZe 



V-u = , 



(2) 



where u is the instantaneous velocity field, P is the 
pressure term, f is an external force per unit of mass, 
the Rossby number is TZo = Uo(2L il)^ 1 (where Uq 
and Lq are respectively normalized velocity and length 
scales taken to be unity) and the Reynolds number is 
TZe = UqLq/v (where v is the kinematic viscosity). 

The forcing term f is introduced in the Navier-Stokes 
equation to study the inverse cascade of energy. In the 



simulations presented in this paper, besides TZe and TZo 
defined at characteristic length scales, we will be inter- 
ested primarily in the Reynolds and Rossby numbers 
based on the forcing scale Lf , at which the external force 
is applied. The latter quantities are defined as follows: 



TZe t 



and 



TZo 



f 



L f U 



U 

2L f n 



(3) 



(4) 



where U is the r.m.s. velocity before the inverse cascade is 
initiated (or equivalently, the r.m.s. velocity at the forc- 
ing scale at any time during the simulation) in units of 
Uq. The time-scale associated with forcing wavenumber 
is defined as, 



T f 



U ' 



(5) 



B. 



Resonant interactions, slow manifold, and 
large-scale structures 



The linear, inviscid approximation of Eq. (nl in the ab- 
sence of forcing has wave solutions called inertial waves 
|10) . In the language of wave turbulence theory, the in- 
viscid version of Eq. ([I]) can be re- written as [24] : 
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where * denotes complex conjugate, Cl'° SpSci — (s q q — 
s p p)(h* s x h* ) ■ h*/2 is the modal transfer coefficient, 
and TZo is assumed to be small and therefore repre- 
sents rapid rotation (hence weak nonlinearity) . Note that 
Eq. (pH) is not closed and, hence, any practical solution 
can be realized based on an equivalent closed set of equa- 
tions (see e.g. [25] for a thorough discussion on closed 
models). Furthermore, the Craya-Herring helical basis 
h s [13 [27] has been used in deriving Eq. ^ with the 
canonical basis corresponding to a given wavevector k, 
■ijj := k x z and k X ■0. The amplitude a s of u is asso- 
ciated with the helical wave with a dispersion relation 
for the wave frequency: w s (fc) = 2fls-^-. Each wave vec- 
tor is associated with two waves of opposite polarization, 
s = ±1. Clearly, tu s (k) = implies a flow restricted to 
a plane perpendicular to the rotation axis (i.e., fcii = 0). 
Hence, 2D modes are also known as slow modes. In other 
words, 2D modes correspond to vortical motions with no 
fast wave modulation. 

The mechanism of transfer of energy towards two di- 
mensional modes that is responsible for the formation of 
Taylor columns is based on near resonant condition of 



the interacting triads 



Sk 



n 
p 



q 



0(Ho) with k + p + q = 0. (7) 



However, the problem with wave turbulence theory is 
that it is not valid for too small values of k\\ . In fact, the 
predicted energy transfer is zero for fen = |28| because 
2D and 3D modes are decoupled in such theories at low- 
est order. Similar analysis is presented using two-point 
closures of turbulence, such as the Eddy Damped Quasi- 
Normal Markovian (EDQNM) closure developed earlier 
in the context of rotating flows (see, e.g., [2H])- Even 
a sophisticated asymptotic quasi-normal Markovian the- 
ory, built on the EDQNM closure [SQlEI], does not deal 
with fen = 0. Thus, while the gradual concentration of 
energy in close proximity of the slow manifold can be the- 
oretically justified to explain numerical and experimental 
observations, the exact coupling between the slow man- 
ifold and the 3D modes leading to a transfer of energy 
from 3D to 2D modes still remains an unresolved prob- 
lem. The inverse cascade of energy, that will be further 
elaborated upon in Sec. |IV[ presumably happens in this 
slow manifold. 

An alternative theory on the egression of columnar 
structures is given by [12]; it is based on the conserva- 
tion of linear momentum P z = H„ (x x uS) z dV and 

of angular momentum L z — J v (x x u) z dV in the ax- 
ial direction (within a cylinder of radius R), resulting 
in a relative concentration of the kinetic energy density 
within this cylinder where it disperses to form columnar 
clouds. This holds in the linear time scale fi" 1 , when the 
non-linear term is small and hence can be neglected in 
comparison with the Coriolis term (Uq <C CILo). How- 
ever, the percentage of total energy contained within the 
cylinder falls as (fit) -1 , so the columns eventually be- 
come weak, although the energy density remains higher 
within the cylinder than outside. The time scale asso- 
ciated with this process, tq ~ fi -1 , will be relevant for 
the analysis of the inverse cascade regime in the following 
sections. 



C. Coupling between modes and energy transfer 

Lewis Fry Richardson's famous couplet, l Big whirls 
have little whirls that feed on their velocity, and little 
whirls have lesser whirls and so on to viscosity' is the an- 
tithesis of observations made by experimentalists [5] and 
analysts [15, 16 in the context of rotating turbulence. 

The notion of inverse cascade of energy to large scales 
is well known in 2D turbulence (32] (also see e.g., [T7] ) 
and may be justified in simple terms on the basis of 
Fj0rtoft's theorem due to the conservation of quadratic 
invariants (see, e.g., [IS]). In other words, nonlinear tri- 
adic interactions conserve both the energy and the en- 
strophy, Z := (oj 2 ) /2, and as the latter is advected to- 
wards smaller scales, a fraction of the energy cascades to- 
wards larger scales to maintain the balance in each triad. 



As discussed in the introduction, the justification for an 
inverse cascade of energy in three-dimensional rotating 
turbulence is not so straightforward since the conserva- 
tion laws change in three dimensions. Nevertheless, a 
similar argument could be made in the case of three di- 
mensional flows based on helicity and the possibility of 
an inverse cascade of energy may be alluded to, as has 
been explained in [531 133] ■ I n fact, it has been argued 
in [21] that interactions between three helical modes of 
the same polarization s, will lead to an inverse energy 
cascade for the same reason as has been postulated by 
Kraichnan [34 and Fj0rtoft in the two dimensional case. 
Numerical simulations of three dimensional flows, where 
the non-linear interactions have been restricted to iden- 
tically polarized wave numbers in all triads, further cor- 
roborate the aforementioned argument |35j . 

In the previous subsection we have summarized theo- 
ries that clearly vindicate the notion of a gradual transfer 
of energy towards the slow manifold, without being able 
to formally account for the exact coupling between the 
3D and the 2D modes. However, once the energy is in 
the 2D modes and if the coupling between the 2D and 
the 3D modes is weak, one can naively expect an in- 
verse cascade to develop as in the case of 2D turbulence. 
The strength of the coupling between the 2D and the 
3D modes has been studied by [36] and also by [37] us- 
ing numerical simulations. In the case of infinitely small 
IZo and in a periodic box, the 2D modes are expected 
to decouple from the 3D modes and evolve under their 
own dynamics. This is in agreement with the evidence of 
decoupling between the 2D and the 3D modes that was 
observed in numerical simulations of freely decaying ro- 
tating turbulence [38] and of ideal helical rotating flows 
[55] , However, note that some authors claim that these 
modes never decouple in infinite domains |40j . 

The decoupling is further illustrated below based on 
the presentation in [37] and extended to consider the 
flux of energy interchanged between the 2D and the 3D 
modes. It is important to note that Refs. [3"TH3"5] studied 
rotating flows in the absence of forcing, thereby making 
a case for analyzing a completely decoupled set of equa- 
tions for the 2D and the 3D modes; however, Ref. [H] 
considers the effect of forcing. 

We write wavenumbers in three-dimensional Fourier 
space using cylindrical coordinates, k = (kj_,kj[), with 
kj. = (k x ,k y ,0) = (p k ,4> k ), k|| = (0,0, k z ) and k = |k|. 
We denote the 2D modes in Fourier space as U2.o(kj_), 
and the 3D or wave modes as u 3 £>(k). Following [3"T] . 
wave vectors are decomposed as follows: 

W k := {k s.t. |k| ^ and k\\^0}, 
V k :={ks.t. |k| ^0and/c,| =0}. 

Then the velocity field u = (u, v, w) can be decomposed 
as: 

n (V\ - / u 3^( k ) iikeW k .. 

uw - 1 u x (k ± ) + u/(kj_)i if k e V k [ *> 



where u 2 £>(kj_) = uj_(kj_) + w(k±)z. Likewise, the to- 
tal energy E — ^ k |u(k)| 2 /2 may be written in terms of 
the modal components as E = E 3D + E 2 d = E 3D + 
(E ± + E w ), where E 3D = J2 ke w k l u ( k )l 2 A E x = 
J2 k± |u 2D (k ± )| 2 /2 and E w = £ k ^ ^(k±)| 2 /2- 

Multiplying the spectral form of Eq. (Til) by u*(k) and 
integrating over all wavenumbers in W k and Vu respec- 
tively, results in two differential equations for the total 
energy in the wave and the slow modes as follows: 

dtE 3 £> = n 2 n^3_D — n 3 £> + e 3 D, (9) 

dtE 2 D = — n 2 D^3D — n 2 £) + C2D, (10) 

where the e 3 o and e 2 ]j terms refer to the corresponding 
components of the forced energy injection, and 3D (2D) 
refers to fen 7^ (kn — 0), as stated before. Equations ([9]) 



and ( 10 1 are congruous to the equations derived in 



When positive, the term H 3 d refers to the 3D energy 
that is transferred to small scales and dissipated per unit 
of time (thus balancing the e 3 £> term) and results from 
triadic interactions that move energy from the 3D modes 
to the 3D modes (resonant interactions involving three 
fast modes or those between two fast modes and one slow 
mode). Similarly, the term H 2 d results from all triadic 
interactions that move energy from the 2D modes to the 
2D modes and, when positive, it's net effect is to balance 
the injection of energy per unit time in the 2D modes. 
Finally, H 2 d^3D is the flux of energy across fen = in 
Fourier space, i.e., energy going from the 2D to the 3D 
modes when H 2 d^ 3 d (t) > 0. This term is expected to 
be 0(Ho) [37], and as a result, in the limit of zero Ho, 
the slow manifold decouples from the wave modes and 
the energy equations ^ and (10 1 are as follows: 



dfE 3 D — — 11323 

dtE 2 o = — n 2 £> 






It may be emphasized that, in this limit, W 2 d only in- 
volves triadic interactions between slow modes and H 3 d 
involves interactions between fast modes. Moreover, the 
equation for the evolution of the 2D energy further de- 
couples into equations for E± and E w (see, e.g., [37]). 



III. NUMERICAL SIMULATIONS AND 
SUBGRID SCALE MODEL 

A. Large eddy simulations 

We integrate the Navier-Stokes Eq. (Il]) in a rotating 
frame of reference using a parallel pseudo-spectral code 
with periodic boundary conditions |42j . A second order 
Runge-Kutta method is used to evolve the equations in 
time and no dealiasing is done because a LES is used; so 
the maximum resolved wavenumber is k c = N/2, where 
N is the linear resolution. As large scale separation be- 
tween the box size and the forcing scale is essential to 
study inverse cascades with reasonably large values of the 



Reynolds number, we use large eddy simulations (LES). 
The sub-grid scale model is such that the wavenumbers 
below a cut-off wavenumber, k c , are resolved explicitly 
whereas larger wavenumbers are modeled based on en- 
ergy and helicity contributions to the eddy viscosity and 
the eddy noise terms in the EDQNM equations. For com- 
pleteness, the model is summarized below. 

First, the larger resolved scales are computed by inte- 
grating the following equation: 



d ,2/1 

1 



u a (k,t) = r<(k,t) 



-^-Pa/3£^ 7 CUc(k,t)+f Q (k,t), (11) 

which is basically the Fourier transform of Eq. (IT]) (bar- 
ring the newly introduced subgrid model term, v k \ ko ). 
Here the greek subindices denote cartesian components 
of the vectors and tensors and Einstein summation con- 
vention is assumed. The term T<(k, t) is the Fourier 
transform of the nonlinear term in Eq. (nl) over all modes 
with k < k c . In other words, it represents the nonlin- 
ear transfer that arises from the convolution sum over 
a truncated triadic domain, k + p + q = 0, k,p,q < k c . 
This term is computed using the pseudo-spectral method. 
The eddy viscosity term, v k \ kc is associated with the sub- 
grid model and is computed based on parameters that are 
modeled from the unresolved scales. P Q , / g(k) = 5 a p — ^- 
is the projector operator on the solenoidal velocity field 
and e is the antisymmetric tensor associated with the 
curl operator (Levi-Civita symbol). 

The isotropic energy spectrum E(k,t) and the helic- 
ity spectrum H{k, t) up to wavenumber 3fc c (including 
unresolved scales) are then obtained through data fit- 
ting and extrapolation from the resolved scales. Next, 
the isotropic energy spectrum E(k, t) for the unresolved 
scales is evolved based on the following: 



(d t + 2vk 2 )E{k,t) 



-2v k \ k k 2 E{k,t)~2v k \ k k 2 H{k,t) 
T P E q (k,t) 



+ T<(k,t) + 



Ank 2 



(12) 



An equivalent balance equation for the unresolved he- 
licity spectrum H(k,t) is solved if the helicity of the 
flow is non-zero (note: when H = , we have v = 0). 
Here, v k \ ka and v k \ kll are terms prescribed by the model 
as before, T^{k,t) represents the energy transferred to 
unresolved scales from the resolved scales and T^ q repre- 
sents the energy and helicity interactions at wavenumbers 
p,q > k c . The analytical forms of the above terms come 
from a two-point analysis of an integro-differential equa- 
tions originating from the EDQNM closure for isotropic 
Navier-Stokes turbulence. Thus, our model assumes that 
isotropy is recovered at sufficiently small scales (smaller 
than the Zeman scale) as was recently shown in a large 
DNS of rotating turbulence [43]. It may be noted here 
that the LES was able to reproduce the results of this 
DNS on a grid of 3072 3 points in which the Zeman scale 



TABLE I: Table of the runs with the total relative helicity 
of the flow pH, the anisotropy exponent j3, the forcing scale 
Rossby and Reynolds numbers, Rof and -Re/, the energy in- 
jection rate e, and the power law index in the inverse cascade 
range of the spectrum of the 2D modes. TG, ABC, RND, and 
ANI respectively stand for Taylor-Green, ABC, random, and 
random anisotropic forcing. Note that pn is the relative he- 
licity of the flow at the time when the inverse cascade starts, 
i.e. at t=0 in the run with rotation. All runs use a grid with 
N = 256 points, a forcing wavenumber kf — 40, an imposed 
rotation Q — 35, and a kinematic viscosity v — 2 x 1CP 4 . 



Run 


PH 


P 


Uof 


Ke f 


c 


index 


TG 


8 x 10 -3 


- 


0.045 


390 


0.030 


»-3 


RND1 


9 x 1(T 3 


- 


0.045 


390 


0.047 


«-3 


RND2 


8x 1(T 2 


- 


0.044 


390 


0.050 


»-3 


RND3 


5x KT 1 


- 


0.046 


420 


0.047 


^-3 


RND4 


7x KT 1 


- 


0.044 


420 


0.047 


»-3 


ANI1 


1 x KT 2 


1 


0.045 


400 


0.010 


£s-3 


ANI2 


8 x 1CT 3 


2 


0.045 


400 


0.010 


»-3 


ANI3 


8 x KT 3 


3 


0.045 


420 


0.007 


« -5/3 


ANI4 


7x KT 1 


3 


0.045 


420 


0.006 


» -5/3 


ABC 


7x KT 1 


- 


0.050 


470 


0.090 


« -5/3 



was resolved [22]. Finally, eddy-noise (upscaling of en- 
ergy towards the resolved scales from the unresolved 
scales) is added to the velocity field based on a recon- 
struction of Eq. (Tl2| . The reader is referred to [20] (see 
Eqs. (27) and (28)) for further explanation. 



B. Description of the runs 

Since the main aim of this paper is to study the in- 
verse cascade of energy, all simulations are forced at high 
wavenumbers kf to ensure sufficient scale separation be- 
tween the forcing and the box scale. We explore different 
forcing functions in order to consider the effects of spec- 
tral anisotropy, number of components in the forcing, and 
the role of helicity on the dynamics of the flow at large 
scales (see table [I]). 

For the Taylor-Green (TG) run, f in Eq. is the 
following [44] : 

fro = /o [sm(k T Gx) cos(k TG y) cos(fc TG z)£ 

- cos(fc TG x) sm(k TG y) cos(k TG z)y] , (13) 

where /o = 5.0 is the forcing amplitude. Such a forc- 
ing function injects zero net helicity in the flow, excites 
only two components of the flow (u and v, although w 
also grows with time as a result of pressure forces) . fxc 
injects energy only into a few 3D modes (no energy is 
injected directly into modes with h\ = 0). The expres- 
sion in Eq. ( 13 ) has many symmetries that are preserved 



during the evolution of the Navier-Stokes equation. To 
break these symmetries and reach a turbulent steady 
state faster, a superposition of two TG forcing functions 
acting at kxa = 21 and 22 was used. Since TG forcing 



involves products of three modes in Fourier space, the ef- 
fective forcing wavenumber is kf = \/3mm{kT G } ~ 40, 
while the projection of the forcing wavenumber into the 
plane of 2D modes is k±j — V2min{kT G } ~ 30. 

In order to study the effect of helicity, we also per- 
formed simulations using the Arn'old-Childress-Beltrami 
(ABC) forcing [45]: 



fABC = fa{[Bcos(k f y) 

+ [Asin(kfx) 4 
+ [Acos(kfx) - 



\- C sm(kfz)] x 
Ccos(kfz)]y 
■Bsm(kfy)]z}, (14) 



with A = 0.9, B = 1 and C = 1.1. ABC forcing is an 
cigenfunction of the curl operator and injects maximum 
helicity (i.e., ¥abc an d V x {abc are co-linear). When 
the flow is forced using this type of forcing function, tur- 
bulence develops only after an instability sets in [1H] • To 
speed up the onset of turbulence, we forced the flow with 
a superposition of two ABC flows at kf = 40 and 41. 
Henceforth, kf refers to the minimum of the two forcing 
wavenumbers. Note that the ABC forcing excites only 
two 2D modes in the Fourier shell with k — kf (in the k x 
and k y axis in Fourier space) and one 3D mode (in the 
fen axis). 

We also used two types of randomly generated isotropic 
forcings. In the first type (labeled RND in table II]), all 
modes in spherical Fourier shells between kf = 40 and 41 
were fired with the same amplitude but random phases. 
The method described in [J7] was used to correlate phases 
and change the helicity of the forcing from zero to maxi- 
mal. This results in isotropic forcing independent of the 
amount of helicity. As a result, more energy is injected 
into 3D modes compared to 2D modes. The second ran- 
dom forcing (labeled ANI in table HI) corresponds to a 
case in which we introduce a new parameter, /3 in order 
to control the extent of anisotropy in the forcing. Wc 
define a new forcing function by multiplying each mode 
in the RND forcing by a factor that concentrates the ef- 
fective forcing near the slow manifold, 



(15) 



f A Ni(k) =[l-f) fnND(k). 



Note that /3 = corresponds to isotropic forcing. Re- 
fer to table [I] for the different values of j3 used in the 
simulations. 

The simulations were started from a flow at rest and 
without rotation and integrated up to ten large-scale 
turnover times (ry = Lf/U ~ 10) until a turbulent 
steady state with a well developed direct energy cascade 
was attained. Next, at a time arbitrarily re-labeled t = 0, 
rotation was turned on and the simulation was continued 
for at least 250 r/ turnover times. A fully developed in- 
verse cascade of energy was observed by this time in all 
the runs mentioned in table |H 



C. Anisotropic spectra 



and 



To study power laws in the resulting inverse cascade 
range of the simulations we refer to both isotropic and 
anisotropic energy and helicity spectra. The decomposi- 
tion of the total energy into energy in 2D and 3D modes, 
Eid and E^d described above can be extended to spec- 
tral densities as follows, based on the definitions pre- 
sented in [35]. 

The isotropic total energy spectrum is computed in the 
simulations as: 



E(k) 



E i"(k)i 2 , 

fc<|k|<fe+l 



(16) 



and is such that the total energy is E = J^k E(k). We 
can also define an axisymmetric energy spectrum: 



e(fcj_,fc||) 



- E 

2 ^ 

fcx<|kxz|<fcj_+l 
fci|<fe 2 <fc||+l 



|u(k)| 2 = e(Mfc), (17) 



where in the latter expression, Ok is the colatitude in 
Fourier space with respect to the rotation axis. The ax- 
isymmetric energy spectrum is such that the total energy 
in 2D modes is E 2 d = J2k ± e (^J-^|| = 0) = J2k e (^> ^k = 
7r/2). As a result, we refer to e(k±, k\\ = 0) as the energy 
spectrum of the 2D modes. 

Reduced perpendicular and parallel spectra can then 
be defined as: 



and 



E(k ± ) = J2< k ^ 



E(k\\) =E e ( fc -L^||) 



(18) 



(19) 



respectively. As for the energy, E — J2k E(kx) — 
J2k E(k\\)- We then introduce the isotropic and per- 
pendicular energy spectra of the 3D modes: 



and 



E 3D (k) = E(k) - e(k, 6 k = tt/2), (20) 



E 3D (k±) = E(k ± ) - e(k ± , fcii - 0). (21) 



Finally, we associate energy fluxes with the energy 
spectra E(k), E(k±), and E(k^). These are defined from 
the transfer functions as follows: 

T(k) = - ]T u*(k) • (tTvu) k , (22) 

fe<|k|<fc+l 



T(k±) = - E u*(k) • (iTw) k , (23) 

k± <|kxz|<fc x+1 



T(fey) = - E u*(k)-(u-Vu) k) (24) 



fcii<fc*<fcii+i 



where the superscript " denotes Fourier transformed 
quantities. Then, the fluxes are as follows: 



n(fc) = - E T(k') 



(25) 



fc'=0 



n (^) = - E T (^)= n ( fc n) = - E T ( fc ii) ( 26 ) 

k' ± =0 fej^O 

These fluxes represent energy per unit of time across 
spheres in Fourier space with radius k, cylinders with ra- 
dius k± and planes with k\\ = constant, respectively. In 
particular, note that II(fc|| = 0) represents energy trans- 
ferred from 2D to 3D modes when positive, and from 3D 
to 2D modes when negative; actually, il2D_>.3_D is the flux 
across slow and fast modes defined in Eqs. ^ and ( 10 ) 
in the previous section. 

Similar definitions can be written for the helicity. 



IV. NUMERICAL RESULTS 

A. Time evolution and spectra 

At the onset, we discuss the results for simulations 
with TG, ABC and RND forcing with and without he- 
licity. The perpendicular spectrum E(k±), the spectrum 
of 3D modes E 3 u(k±) and the spectrum of 2D modes 
e(fcj_, fcj| — 0) at late times are all shown in Fig.lTlfor the 
TG, ABC, RND1 and RND4 runs. The ABC run shows 
a spectrum e(k±,k^ = 0) ~ fc ± and E(k±) ~ fc _1 
large scales. All the other runs have e(k±, h\ = 0) ~ k 



at 

-3 



The time evolution of the spectrum of the 2D energy 
e(k±, fen = 0) for the same runs as in Fig. fllis shown in 
Fig. [21 from t = to 250r/ turn-over times at intervals 
of roughly 35. 7t/ turn-over times. 

The build up of energy at large scales observed in the 
spectra is associated with an inverse cascade of 2D energy 
in the presence of rotation. This can be verified from the 
energy flux that shows a positive range at wavenumbers 
larger than kj (associated with a direct cascade of en- 
ergy) and a negative range at wavenumbers smaller than 
kf (associated with the inverse cascade). Figure pj] shows 
n(jfex) for runs TG, ABC, and RND1. The same behav- 
ior is observed in the isotropic flux H(k) (not shown). 

The fluxes in Fig. [3] do not discriminate between 2D 
and 3D modes, so although they confirm an inverse en- 
ergy transfer, they are not enough to identify what modes 
are responsible for the large-scale pile up of energy. Fig- 
ure [4] shows the time evolution of E, E2D, and E 3 £> for 
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FIG. 1: E(k±), E[io(k±) and e(fcx,fc|| = 0) at late times for 
TG, ABC, RND1 and RND4 forcing (from top to bottom). 



FIG. 2: Time evolution of the spectrum of the 2D energy 
e(fcx,ft|| = 0) for run TG, ABC, RND1 and RND4 (top to 
bottom) upto 250 turn-over times at intervals of 35.7 turn- 
over times. 
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FIG. 3: Energy flux n(Jfcx) for runs TG, ABC and RND1 
(from top to bottom). Solid lines are time averaged while 
dashed are instantaneous fluxes at late times. The fluxes are 
normalized to the value at the forcing wavenumber. 



several runs. In all cases, E 2 d grows monotonically in 
time after a short transient, thereby driving a growth of 
the total energy E. However, the energy in 3D modes, 
E 3 d can either slowly increase or decrease, depending 
on the forcing. On the other hand, the enstrophy in 2D 
modes: 



Z 



2D 



kev fc 



(27) 



remains approximately constant once the inverse cascade 



FIG. 4: Time evolution of E, E 2D , and E 3D in runs TG, ABC 
and RND1 (from top to bottom). 



starts, implying the small scales have reached a steady 
state. The enstrophy in 3D modes, 



Z3D = 



1 



Y^ 



2 ^ 

k£W k 



(28) 



as well as the total enstrophy Z, increases or decreases de- 
pending on the forcing function. Note that the transient 
regime at early times corresponds to a wave dominated 
regime, and is longer for TG forcing as already reported 
in 0S]. Indeed, in the TG case, oscillations in E^£>(t) 
and Zsuit) with the frequency of inertial waves can be 
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FIG. 5: Time evolution of the total enstrophy Z, enstrophy 
in 2D modes Zzd, and enstrophy in 3D modes Z^d in runs 
TG, ABC and RND1 (from top to bottom). 



FIG. 6: Energy spectra at late times for runs ANI1 (top), 
ANI3 (middle) and ANI4 (bottom). 



seen clearly before t/rj ~ 20. 



B. The effect of anisotropic energy injection 

Except for one case (ABC forcing), all simulations in 
the previous subsection seem to show an inverse cascade 
of 2D energy with a kj 3 scaling. What is the origin 

— 5/3 

of the KKBL-likc ~ k ± spectrum in the ABC run? 
Previous studies obtained KKBL scaling with elongated 
boxes jin| or when all triadic interactions between 2D and 
3D modes were shut down (18) (which in fact corresponds 



to the case of KKBL phenomenology). However, in our 
case we used a box with fixed unit aspect ratio and with 
all triadic interactions and coupling between modes were 
accounted for in the simulations. 

Helicity is the most conspicuous property of the flow 
generated by ABC forcing and it is known to affect the 
direct cascade range in rotating flows [51 2H]. However, it 
is predominantly transferred to smaller scales as shown 
later (see also [50J) and runs with isotropic but helical 
forcing (RND4) also show ~ fcj 3 instead of KKBL scal- 
ing. 

The difference between ABC forcing and the other forc- 
ing functions is that fABC excites two 2D modes and one 
3D mode in Fourier space, therefore effectively injecting 
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FIG. 7: Time evolution of the energy spectrum of 2D modes 
e(fcx,fc|| = 0) (top) and of the perpendicular spectrum E(k±) 
(bottom) in run ANI3. A —5/3 slope is indicated as a refer- 



FIG. 8: Helicity spectra H(k±) and h(k±,k\\ = 0) (i.e., the 
spectrum of helicity contained in purely horizontal motions) 
for runs ABC (above) and ANI4 (below) . Slopes are indicated 
as a reference. 



more energy in 2D modes than in 3D modes. The other 
forcing functions considered in the previous subsection 
inject either more energy in 3D modes (RND), or energy 
only in a few 3D modes with no 2D injection (TG). Here 
we explore if anisotropic injection can be responsible for 
the different scaling laws observed in the e(k±,kn = 0) 
spectrum by means of numerical simulations in which we 
control the anisotropy of the external forcing (runs ANI 
in table jl}. 

It should be pointed out that what we call here 
"anisotropic injection" for sake of brevity, is actually a 
more subtle effect associated with how much energy is di- 
rectly injected in 2D modes compared to that into the 3D 
modes. Indeed, TG forcing is directionally anisotropic (in 
the sense that only a few modes in a spherical shell are 
excited in Fourier space) , but it shows similar inverse cas- 
cade scaling as the RND runs, which have no directional 
anisotropy. 

The energy spectra e(/cj_, fc|| = 0), E(k±), and E 3 £>(k±) 
at late times for runs ANI1, ANI3, and ANI4 are shown in 
Fig. [6] The run ANI1 corresponds to a run with random 
forcing with anisotropic exponent ,9 = 1. Runs ANI3 and 
ANI4 have /3 = 3, with zero and close to maximal helicity 
injection, respectively (see table II]). The spectra of ANI2 
behave as ANI1 and are not shown. 



The time evolution of e(k±, kn — 0) and E(k±) for run 
ANI3, up to 250r/ turn-over times at intervals of 35. 7r/ 
turn-over times, is shown in Fig. W\ A clear build up of 
a ~ k ± spectrum can be observed in the runs with 
anisotropic forcing irrespective of the amount of helicity 
in the flow. 



C. Helicity at large scales 

Fig.[8|shows the helicity spectra H(k±) and h(k±, k\\ = 
0) for trie runs with helical forcing ABC and ANI4. There 
is no significant large scale growth of helicity and the rel- 
ative helicity remains negligibly small at low wavenum- 
bers, p H (k) = H{k)/kE{k) -> for k -> (also 
PH(k±) — > for k± — > 0). This is consistent with the 
observation above that helicity does not seem to affect 
the energy scaling in the inverse cascade range, as has 
also been shown in previous studies of helicity cascading 
to smaller scales in rotating flows [5U] . 
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FIG. 9: n(fc||) for runs TG, ABC, ANI1, ANI2 and ANI3 (top 
to bottom). 



D. Coupling and fluxes between slow and fast 
modes 



How does the amount of energy injected into the 2D 
and the 3D modes affect the inverse cascade? Fig. [9] 
shows the flux n(fc||) for runs TG, ABC, ANI1, ANI2, 
and ANI3 specified in table [I] In most runs, the flux is 



negative for small values of fcii (indicating that energy 
goes from the 3D modes towards the 2D modes for larger 
scales), and positive for large values of fen (indicating that 
energy goes away from the 2D modes for smaller scales) . 
Note that as more energy is injected into the 2D modes 
(e.g., as /3 is increased in the ANI runs), the wavenumber 
at which the fluxes change sign moves towards fen = 0, 
and for the ABC flow the flux n(fc|j) is positive every- 
where (i.e., energy goes from the 2D modes to the 3D 
modes at all scales). As a reference, a schematic repre- 
sentation of the excited modes in axisymmetric Fourier 
space (fej_, fen) is shown in Fig. 



10 



The above observations imply that the ABC flow cor- 
responds to the limiting case in which most of the energy 
is injected into the 2D modes and as a result of the im- 
balance, an excess of energy "leaks" from the 2D modes 
to the 3D modes. This can be verified by computing each 
term i n th e energy balance Eqs. Ety and ( 10 ) (see tables 
[TT| and III ) . The flux of 3D and 2D energy can be esti- 
mated from Eqs. ^ and ( 10 ) as follows. Since all terms 
in these equations, with the exception of II2_d and TI3D, 
are known, the equations can be re-written as: 



and 



Tjl.h.S 



n(fc|| =0,t)+e 3D ~d t E 3D , 



nfcs" = -n( fc n = o,t) + e 2D - d t E 2D . 



(29) 



(30) 



The superscript "l.h.s." here and in the table indicates 
that the fluxes are obtained by solving for the l.h.s. of 
the balance equations above. Another way of estimating 
II3C is based on geometrical consideration of the fluxes 
in spectral space (see Fig. 10), 



IT 



:\D 



\max kj < k < kf+1 {U{k,t)} iff = fABC 
]maxy k {U(k) -n(fcjj} if f = fTGjfRND- 

(31) 
Finally, an alternative interpretation of II 3 £> is that it 
transfers energy to the 3D modes with larger wavenum- 
bers where it is eventually dissipated, and hence it bal- 
ances the injection term with dissipation of energy per 
unit of time. Hence, it can also be approximated by 
2v L ma!C k 2 E 3 o{k)dk. The two estimates have been 
found to be of the same order of magnitude in all the 
runs as shown in table |TTJ The fact that the three esti- 
mates are positive and of the same order indicate that 
energy in the 3D modes does not only go to the slow 
manifold but goes to smaller scales where it dissipates. 
The fact that n 2 £> is negative in the ABC run (see ta- 
ble III) is evidence of an inverse cascade of energy in 



the slow manifold once the energy from the 3D modes is 
transferred to the 2D modes. For the other runs, even 
though n 2 £> is positive and small in magnitude, it merely 
implies that more energy cascades to the smaller scales 
than to the larger scales. It is important to note that 
the positiveness of H2D hints at positive eddy viscosity 
and the possibility of the inverse cascade of energy in the 
slow manifold cannot be ruled out. In fact, the nature of 
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TABLE II: Amplitude of the terms in Eq. pM. The time derivative dE^ujdt was obtained using centered finite differences from 
the data. II(fc|| = 0) represents energy per unit of time transferred from 2D to 3D modes, and e^D is the power injected in 
the 3D modes, Illjjj 8 ' is the flux of energy in the 3D modes estimated from Eq. (29 1, II^' 
2v J k 2 Z3o{k) dk is an estimation based on the energy dissipation rate. 



is estimated from Eq. (311, and 



Run 



dE- iD /dt 
1.0 x 10" 4 
1.0 x 10~ 4 
4.0 x 10~ 5 
1.0 x 10" 4 
1.0 x 10~ 4 
1.0 x 10~ 4 
3.0 x 10~ 5 
-2.0 x 10~ 4 



n(fe|| = 0) 



£313 



TT- 
11 j 



2iyJk 2 E 3 p{k)dk 
1.0 x 10~ 2 
2.0 x 10~ 2 
3.0 x 10~ 2 
4.0 x 10" 2 
1.0 x 10" 2 
5.0 x 10" 3 
4.0 x 10" 3 
3.0 x 10" 2 



TG 
RND1 
RND4 
ANI1 
ANI2 
ANI3 
ANI4 
ABC 



-2.0 x 10" 3 


3.0 x 10~ 2 


2.8 x 10~ 2 


1.0 x 10 


-6.8 x 10~ 3 


4.6 x 10~ 2 


3.9 x 10" 2 


2.0 x 10 


-6.3 x 10" 3 


4.6 x 10~ 2 


3.9 x 10~ 2 


2.0 x 10 


-5.2 x 10~ 3 


8.9 x 10" 3 


3.6 x 10" 3 


4.0 x 10 


-2.3 x 10~ 3 


9.3 x 10~ 3 


6.9 x 10~ 3 


2.0 x 10 


-6.9 x lO" 4 


6.5 x 10" 3 


5.7 x 10~ 3 


2.0 x 10 


-6.0 x 10~ 4 


5.4 x 10~ 3 


4.7 x 10~ 3 


2.0 x 10 


8.3 x 10~ 2 


2.0 x 10~ 2 


1.0 x 10" 1 


3.0 x 10 



TABLE III: Amplitude of the terms in Eq. (10 1. The time 



derivative dE2p/dt was obtained using centered finite differ- 
ences, €2D is the power injected in the 2D modes, and II^' 8 ' 
is the flux of energy in 2D modes estimated from Eq. (301. 



Run 



TG 
RND1 
RND4 
ANI1 
ANI2 
ANI3 
ANI4 
ABC 



dE 2D /dt 
4.0 x 10" 4 
1.0 x 10" 3 
2.0 x 10" 3 
1.0 x 10" 3 
4.0 x 10~ 4 
7.0 x 10~ 5 
8.0 x 10~ 5 
5.0 x 10~ 4 



£2D 



1.0 x 10" iu 

1.3 x 10~ 3 
1.5 x 10" 3 

1.1 x 10~ 3 
7.0 x 10" 4 
5.0 x 10~ 4 

6.4 x 10" 4 
7.0 x 10~ 2 
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10" 
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10 _ 
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10" 
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the evolution of the e(kj_,k\\ — 0) spectra over time, as 
shown in Fig. [21 is evidence of occurrence of the inverse 
cascade of energy in the slow manifold. It may also be 
worth pointing out that with increasing j3 (anisotropy), 
IT2D, for the ANI runs, become less positive and seems 
to approach the nature of the energy cascade exhibited 
by the ABC run (see table |JTij). 

The picture that emerges for the fluxes from the values 
in tables [LT] and |HI| is illustrated schematically in Fig. [TO] 
For isotropic forcing, a fraction of the energy injected into 
the 3D modes is transferred to the slow manifold and the 
remaining 3D energy is transferred to the 3D modes with 
larger wavenumbers. 

Recently, |H] did a detailed study of energy transfers 
in forced rotating turbulence and concluded that the for- 
mer transfer, from the 3D to the 2D modes, is non-local. 
The results in [JT] (with simulations forced at smaller 
wavenumbers than in our case) are consistent with our 
results, except that they attribute the development of the 
~ fc~ 3 spectrum in the 2D modes to a direct cascade of 
enstrophy once the 3D modes inject energy directly into 
the 2D modes with the smallest wavenumbers. Although 
our analysis cannot disprove this conjecture, the ampli- 
tude of the terms in tables HI] and IIHI and the nature of 
the evolution of the 2D energy spectra in Fig. [2] hints at 
a likely inverse transfer of energy in the slow manifold. 



FIG. 10: A schematic depiction of the direction of transfer of 
energy (and corresponding fluxes) in the case of forced rotat- 
ing turbulence. Here f m is the normalized unit amplitude of 
forcing in Fourier space. The black dots indicate the modes 
that are directly excited by the different forcing functions, 
with 2/ m indicating twice the energy injected in that mode 
compared to the energy injected into other modes. Directions 
of the arrows are based on the analysis of data from tables III] 
[Inland Figs. [2] [4] 



For anisotropic and ABC forcing the picture changes. 
As more energy is injected directly into the 2D modes 
by the forcing, the flux of energy from the 2D to the 
3D modes, H(At|| = 0) in table ITT) increases (from larger 
negative values to smaller negative values) and eventually 
reverses sign becoming positive. In the ABC flow (and 
likely for other flows with very high anisotropic forcing) , 
the energy injected directly into the 2D modes undergoes 
an inverse cascade in the slow manifold, and later the 
excess of energy in these modes relative to the 3D modes 
"leaks" energy into the 3D modes at large scales (see 
Fig. FLO]). As the effect of the 3D modes over the 2D 
modes is less relevant in these runs, the runs with either 



small or positive n(feii = 0) display 



-5/3 



scaling. 



It is also interesting to point out that the 2D and 
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the 3D modes are not necessarily decoupled, as can be 
seen, e.g., by comparing the ratios n(/c|| = 0)/U^j s - and 
II(fc|| = O)/!!!;^)' 8 '. The relevant quantity for either a 



k 



-3 



-5/3 



scaling in the inverse cascade energy 
spectrum is n(feii = 0). Indeed, based on the previous 
discussion, if little energy goes into the 2D modes from 
the 3D modes, or if energy goes from the 2D modes into 
the 3D modes with the 2D modes being the most ener- 
getic and dominating the dynamics, we can assume that 
the cascade in the slow manifold is dominated by the 
turnover time tj_ ~ lj_/uj_ (where l± is a characteristic 
length in the slow manifold and u± the 2D r.m.s. velocity 
at that length). With only one relevant timescale, KKBL 
phenomenology tells us that the energy flux goes as: 



n 



2D 



T± 



(32) 



which results in a 



-5/3 



scaling law. On the other 



hand, if energy goes from the 3D modes to the 2D modes, 
interactions with the 3D modes cannot be neglected. Be- 
sides the slow turnover time t± , we have to consider now 
the 3D turnover time and the timescale associated with 
the fast waves, tq ~ ^. There is no unique dimensional 
solution in this case but we can borrow from the phe- 
nomenology developed by Kraichnan |51j for magneto- 
hydrodynamic (MHD) turbulence where the effect of the 
waves modulates the dominant time-scale of the flow. 
This phenomenology has been successfully extended to 
rotating flows 52, 53 , including in the helical case [6]. It 
states that in the presence of waves, the nonlinear trans- 
fer is slowed down because of the waves and the relevant 
parameter of the problem is the Rossby number, i.e. the 
ratio of time-scales of the wave and the nonlinear turn- 
over time. Thus, we can assume that the flux will be 
slowed down by a factor proportional to the ratio tq/t 
where r is the relevant (and unknown) turnover time for 
the problem, 



n 



u 



2D 



T T 



(33) 



It is interesting to point out that if the turnover time 
in the above expression is built upon the velocity at the 
forcing scale f/j as r ~ lj_/V / (i.e., assuming interac- 
tions are non-local in the inverse cascade range and that 
most of the energy in the slow manifold comes directly 
from the 3D forced modes, which is consistent with the 
large and negative values of II(fc|| = 0) in some of the 
runs), Eq. ( [33] ) results ina~ kj 3 scaling for the energy 
spectrum of the 2D modes. Note that such a choice of 
the time-scale is consistent with the non-local transfers 
reported in |41j , and that non locality of nonlinear trans- 
fer in rotating turbulence in the direct cascade was also 
observed in 1351. 
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FIG. 11: Spectrum of the maximum eigenvalue of the rate 
of strain tensor in horizontal planes, for different runs as a 
function of time. The amplitudes of the spectra are normal- 
ized to their initial values (when the rotation is turned on 
corresponding to 77 = 8). From top to bottom: TG, RND1, 
RND4, and ABC. 
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E. Large-scale shear 



CONCLUDING REMARKS 



The transfer of large-scale energy from the 2D modes 
to the 3D modes in the ABC run should have an impact 
in the large scale structures that develop as a result of the 
inverse cascade. In this section, we show that large-scales 
in the ABC run have large shear, while in the other runs 
(in which energy is mostly transferred from the 3D to the 
2D modes or cases where the transfer is negligible) have 
much smaller shear at large scales (see Fig. [TT| ). The 
development of large-scale shear in the ABC run also 
introduces a new timescale for the 3D modes, r s /j which 
is independent of the scale and therefore consistent with 
the observed ~ fc _1 scaling in the total energy spectrum. 

To study the effect of shear, we analyze the velocity 
gradient tensor, 03 defined as follows: 



2J:= Vu 



d x u d y u d z u 
d x v d y v d z v 
d x w d y w d z w. 



(34) 



The velocity gradient tensor may be written as the 
sum of the symmetric rate of strain tensor, 6 and the 
antisymmetric rotation tensor, 9^, i.e. 2J = & + 9t. Note 
that & may also be written as: 



6 = ^(» + 23 T ) 



(35) 



We analyze the spectrum of the maximum eigenvalue 
Araaz of the rate of strain tensor 6 in a horizontal plane 
for several runs (in Fig. 11 note the amplitude of the 



spectra of \ m ax in the figure is normalized by its initial 
value before rotation is turned on). In most of the runs, 
shear decreases when rotation is turned on and the spec- 
trum seems to reach a steady state (with significantly 
less shear than the isotropic and homogeneous turbulent 
flow) at late times. However, in the ABC run, shear in- 
creases at large scales as time evolves (with a decrease in 
shear at small scales). The increase in large scale shear 
in this run can be understood in the light of the previous 
discussion and based on the fact that n(fcii = 0) > 0. 
In the case of ABC forcing, 2/3 of the energy injection 
corresponds to the 2D modes and the energy in the slow 
manifold undergoes an inverse cascade. However, some 
of the large-scale energy in the 2D modes is transferred 
back to the 3D modes (as is evident from n(fcii = 0) > 0) 
and the excitation of large-scale 3D modes in turn creates 
large-scale shear (note: d z {-) = in the slow manifold 
when fcii = but d z (-) ^ in the 3D modes). 

It is interesting that once large-scale shear is present, 
a new constant time scale (i.e., independent of length 
scale) appears. Large scale shear is associated with the 



shear time scale, t s 



where X max (x,y, z) is 



max{A maI } ' 

the largest eigenvalue (in magnitude) of 6 at any given 
location (x,y,z). Dimensional analysis then hints at a 
flatter energy spectrum, and a fc _1 energy spectrum (as 
observed for the total energy in the ABC run) is not 
uncommon in shear dominated flows [54 . 



An attempt has been made to provide a coherent treat- 
ment of the various large scale physical processes involved 
in rotating flows with a special emphasis on the break- 
ing of universality in the inverse cascade of energy for 
rotating turbulence due to the anisotropy in the forcing. 

We observe a ~ fcj spectrum in the 2D modes when 
the forcing is isotropic and the observed spectrum can 
be associated with an inverse transfer of energy. This 
spectrum has been reported before [T4l BTj although it 
must be said that in [41], the observed kj 3 spectrum is 
attributed to a direct enstrophy cascade, as is thought 
to happen in the case of two-dimensional Navier-Stokes 
turbulence (see (25] and references therein). 

We also find that in the case of strongly anisotropic 
forcing such as the ABC, the inverse cascade for the 2D 

— 5/3 

energy follows a k ± ' spectral law but the 3D energy 
follows a shallower law, E(k) ~ k~ x attributed to the 
creation of shear at large scales. Of course, a A: -1 spec- 
trum is not unheard of in turbulence; it arises for a field 
which is advected at a constant rate, e.g., for a passive 
tracer in a turbulent flow and it is also documented in 
shear flows [55] (see also the recent review in [56]). Al- 
though a recent analysis of atmospheric data indicates 
that the energy cascade is forward at all scales (see [57] 
and references therein), it should be noted that the fc -3 
and k~ 5 ' 3 spectra are both observed. 

There is no indication of inverse cascade for the helic- 
ity in any of the runs. However, it is not clear what the 
origin is of the rather flatter spectrum of the helicity at 
large scales (see Fig. [8] and also [22], Fig. 3). It is con- 
sistent with the fact that the relative helicity becomes 
negligible at large scales, so this could be simply inter- 
preted as eddy-noise. This point will have to be investi- 
gated further. Note that a k~ l scaling law, which implies 
less energy at large scales, may be consistent with toned 
down 3D non-linear dynamics due to the strong helicity 
in the flow. 

Finally, it must be said that these results are obtained 
with a sub-grid scale model and due to the non-local in- 
teractions that seem to be present in the inverse cascade, 
the modeling of the small-scale may affect the behavior of 
the large scales. This effect is perhaps of less consequence 
when compared with models using an hyper-viscous term 
with the dissipation proportional to ~ k 2a E(k). Indeed, 
in our LES the computations of the eddy viscosity and 
eddy noise take into account the energy and the helicity 
spectra up to 3fe c , where k c is the cut-off wavenumber. 
We would like to mention that DNS of the flows presented 
in this paper are being performed currently to substan- 
tiate the observations based on the LES. 
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